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Singlet-triplet spin qubits in six-electron double quantum dots, in moderate magnetic fields, can 
show superior immunity to charge noise. This immunity results from the symmetry of orbitals in the 
second energy shell of circular quantum dots: singlet and triplet states in this shell have identical 
change distributions. Our phase-gate simulations, which include 1// charge noise from fluctuating 
traps, show that this symmetry is most effectively exploited if the gate operation switches rapidly 
between sweet spots deep in the (3,3) and (4,2) charge stability regions; fidelities very close to one 
are predicted if sub-nanosecond switching can be performed. 



The spin degree of freedom of the few-electron quan- 
tum dot (QD) is an excellent building block for a qubit. 
While a single electron spin may serve directly as a qubit 
[T] , the difficulty of single-qubit operations makes it de- 
sirable to encode a qubit in a multi-electron state. Con- 
siderable success has been achieved with a two-electron 
encoding [2] , in which the singlet and spinless triplet lev- 
els of the double quantum dot (DQD) define a logical 
qubit [3J. Electric pulses, applied on the microsecond 
scale, permit all necessary one- |H \5\ and two-qubit [1T| E] 
operations when supplemented by magnetic field gradi- 
ents [MTU]. 

This paper addresses problems with the crucial "ex- 
change gate", which, while it has provided a route to im- 
pressive progress in the singlet-triplet qubit, suffers con- 
siderably from low-frequency noise due to charge traps 
PTTHT5] . In this gate a DQD is moved away from the 
exactly "neutral" electron distribution, i.e. having one 
electron on each QD (referred to as (1, 1)) to one having 
a slight bias towards double occupancy of one QD (e.g. 
the left QD: (2,0)). Only the singlet configuration per- 
mits electron transfer from (1, 1) to (2, 0), while transfer 
from the triplet state is blocked (Pauli spin blockade). 

Here we show that, paradoxically, the exchange gate 
will be much less susceptible to change noise if the DQD 
is pulsed fully from the (1, 1) to the (2, 0) regime. Pulsing 
far into the (2, 0) region lifts the spin blockade for the 
triplet as well as the singlet state [HI , making the singlet- 
triplet splitting highly protected from charge noise. We 
show that the fidelity of this gate will be excellent under 
two conditions: 1) the pulse rise and fall times should 
be sub-nanosecond; and 2) the electrons should be in 
the second shell, so that singlet and triplet states have 
the same charge distribution. This means that the best 
exchange gate is predicted to occur for the six-electron 
DQD with four non-participating "core" electrons, so that 
the desired transition is actually between (3, 3) and (4, 2). 

Model — Our description of DQDs starts with the 
single-particle eigenstates of a circular QD with confining 
potential V (x, y) = '^r 2 [15JDI] and out of plane mag- 
netic field B. The eigenstates are the Fock-Darwin (FD) 



states ip at p [T7], with energies E a> p — (2a + \f3\ + 1) Ml— 
Phw c HUES]; w c = and tt 2 = uj 2 + u? c . We will con- 
sider moderate B fields: the degeneracies E a< p, for the 
same /?, are lifted; but E ai p of different a do not cross 
(ujc/uo <C 1). The single-particle eigenstates are grouped 
into "atomic" energy shells [20] • The ground state ipo,o is 
well separated from the first two excited states "00, ±i- 

We employ a description for few electron DQDs that 
takes into account multiple energy levels and electron- 
electron interactions [211 [22] . As in the work of Burkard 
et al. [15], we use the FD single-particle eigenstates 
as the basis states of a Hubbard model. For two- 
electron DQDs, we include only the (1,1), (2,0) and 
(0, 2) electron configurations. The singlet (S) and s z = 
triplet (T) state can be written as the product of its 
spin and orbital part: * S/T = {<j> u (j) 2 } s/a ® ln) ^ T> . 
The electrons occupy states (pi, which need to be sym- 
metrized/antisymmetrized for the S/T-state (as indi- 
cated by {»,»} s / a ). In the (1,1) configurations, <f>x/4>2 
is close to the FD ground state ip ' on the left /right 
QD. In the (2, 0) and (0, 2) singlet configurations, both 
electrons fill the same orbital ground state, close to iPqq R 
on the respective QD. For the triplet, the Pauli exclusion 
principle requires that the next higher orbital state ipo { R 
is also occupied. 

As in atoms, the first electron shell ipob R is completed 
with two electrons. We assume that in the six-electron 
configuration the first two electrons on each QD are 
paired in a singlet state. We thus adopt a frozen core 
approximation: The (3, 3) configuration for six-electron 
DQDs is therefore represented as a (1, 1) valence config- 
uration (and similarly the (4, 2)/(2,4) and (2,0)/(0,2) 
configurations). One just needs to use the appropriate 
orbital wave function of these "valence" electrons. The 
valence orbital ground state is then V'oi^i while the first 

excited state is 4>o^-l- 

The DQD Hamiltonian is expressed in the basis 
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The diagonal entries describe the energy of each state. 
The difference between (1,1) s and (1,1) T matrix ele- 
ments is neglected, since it is commonly small [E]. Un- 
equally occupied QDs are higher in energy by Us/t |23| . 
A = Us — Ut is the energy difference between the doubly 
occupied states. Electrostatic bias, modeled by the pa- 
rameter e, influences the relative state energies of uniform 
and unequal electron arrangements. The off-diagonal el- 
ements in Eq. (JT|) describe the spin-conserving hopping 
process of electrons between the dots. 

Fig. [T] shows the energy spectrum as a function of e. 
Close to state degeneracies |e| = U a , the hopping process 
hybridizes electron configurations of the same total spin. 
The ground state Es/Et is shown in blue/red. At e = 0, 
both energy levels are mainly in the (1,1) charge config- 
uration, and their energy difference is minimal. Es and 
Et are lowered in energy for increasing bias, due to the 
transfer of electrons between the QDs. For large e, the 
ground states are close to (2,0) ST with an energy differ- 
ence A; we indicate one point deep in the (2, 0) region as 
the "high-bias" configuration e = €hb- 

Our treatment of few-electron DQD is not self- 
consistent; it employs energy spectra of single-particle 
states, which are successively filled with electrons. The 
FD-states are a valid ansatz for the description of few 
electron QDs, if the electron-electron interaction influ- 
ences the single-particle energies weakly or shifts all en- 
ergy levels by a fixed value. The last scenario is consis- 
tent with the calculations of Guclu et al. , where the addi- 
tion energy of interacting electrons has a constant offset 
compared to the non-interacting case [25]. This is con- 
sistent with the Constant Interaction Model, introduced 
by Averin and Likharev |261 127] , in which the energy 
spectrum of QDs remains unchanged when an electron 
is added to or removed from a QD. 

Charge Noise — Charge noise is generally modeled by 
a random distribution of classical two level fluctuators 
(TLF), which couple electrostatically to QDs HE1I2H]- If 
the occupations of the charge traps (CT) vary with a 
broad distribution of fluctuation rates, 1/f noise is gen- 
erated. The coherence of the QD decreases, as seen by 
the time evolution of superpositions: 
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Est (t) is the time- varying energy difference of the qubit 
levels, which deviates from the ideal value due to the 
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Figure 1. Energy diagram for two- and six-electron DQDs, 
as described by Eq. Q: Electrostatic bias, modeled by e, 
transfers one electron from uniform electron distribution on 
the two QDs towards two excess electrons on the left QD. 
The blue/red line represents the singlet/triplet ground state 
Eg/T', black curves are excited states. Charge noise generates 
fluctuations between Es and Et, as described in the main 
text. The electron configurations are highly insensitive to 
charge noise at e = (the "neutral" electron configuration) 
and e — eHB (the "high-bias" configuration far away from 
the two anticrossings) . The inset shows the charge stability 
diagram following van der Wiel et al. [24J. (n_t,n_R) is the 
stable charge configuration on the left and the right QD. Vl/r 
describes electrostatic voltages applied to the left /right QD; 
red arrows indicate gate tunings corresponding to the energy 
diagram. 



coupling to TLFs: SEst = Est — (Est)- (• • • ) describes 
averaging over many experiments. Assuming a static en- 
vironment during one run, the coherence time Ti is re- 
lated to the statistics of the TLFs: 2TT 1 = er^ /V^r; 

°se = Coo dtsE 2 (t)m mm . 

We analyze the relative energy shift of the qubit levels 
of a DQD which couples to a CT. In first order perturba- 
tion theory the fluctuation of the singlet-triplet splitting 
is described by [33 : 



5E& = 



* T \e$>cT\ *t> - (*S \e<S>cr\ *s> 



(3) 



$ct is the electrostatic potential of a CT. Since for QDs 
that are suitable for qubits, CTs are at some distance 
from the QD center, we make a multipole expansion of 
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ro is the position of the CT relative to the center of the 
DQD, r is the QD electron coordinate. This expansion 
resolves the coupling of a TLF into dipole (— E • d) and 
quadrupole (-(1/2) {d l E J )-Q^) terms; <f = e(* \r l \ 
and Q l: > — e \r l r 3 \ ^) are the first two electric mo- 
ments of the DQD. We analyze two points in the charge 
stability diagram e = and e = 6hb ("sweet spots", intro- 
duced in Fig. [T]) at which coupling is weak to TLFs. High 
couplings are obtained if the qubit states have different 
dipole moments, which generate energy shifts scaling like, 
e.g., l/rl. 
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The eigenstates of the singlet-triplet qubit of Eq. ([I]) 
can be approximated at e = 0: \S/T) cx |11) S / T — 

IfJ^ (\ 0, 2) S / T + | 2, 0) g/rj ■ These have equivalent 
dipole moments for the two qubit levels; the charge dis- 
tribution of a DQD arranged in x-direction has mirror 
symmetry to the y-z plane. The quadrupole contribu- 
tion describes the spread of the charge distributions. 
The unequal degree of hybridization of the singlet and 
the triplet state creates different variances in x-direction: 

- (0) 2 - fe) 2 ) • e -f • ^ Of - 4) • The 

first factor describes the hybridizations for the singlet 
and the triplet, the second factor involves to the inter- 
dot distance do of the DQD, and the third factor is the 
gradient contribution of the electric field of the TLF. It 
describes an energy shift proportional to the hybridiza- 
tion of the ground state {t s /t /Us/t) > which decays like 
1/rg in the TLF-QD distance. A similar expression holds 
for the six-electron DQD. 

Considering the two-electron DQD for high bias (e = 
chb), the left QD is lower in energy than the right QD. 
Since the charge configuration is (2,0), the qubit states 
are close to the FD-states {V'ooj TOi} ■ The dipole contri- 
butions to Eq. (|3| vanish. The quadrupole contribution 

of Eq. § is: SE^ « (^) ; ^ - 2|). The 

first factor comes from the different spread of the density 
of the qubit states, while the second factor describes the 
influence of the CT. We find a 1/rg-scaling in the CT-QD 
distance as for the low-bias sweet spot. 

The situation improves for six-electron DQDs. As 
the valence electrons' wave functions Vfyii 316 com- 
plex conjugates of each other, not only the quadrupole 
term of Eq. ([3|, but all multipole contributions van- 
ish. SEgj, depends on the charge density of the single- 
electron wave functions, as e&cT describes interactions 
with only one of the two valence electrons. The second- 
order dipole contribution of TLFs (second-order Stark 
effect) vanishes accordingly, since it involves only an 
overall shift of the confining potential. The first non- 
vanishing contributions are second-order quadruple cou- 

plings: 8E§ « -§ {d iEj f ^gfe^ - We note 
that this contribution has 1/rg scaling with the CT- 

(2) 

QD distance, which suppresses SE ST considerably. This 
protection criterion for six-electron DQD is strongest 
for perfect circular symmetry. For weakly elliptic QDs, 

2 

V = — p (X -\- p cos (2</>) ) , the diagonal terms of the 
quadrupole tensor differ, weighted by the ellipticity f3: 

Q „ _ QTg ,yy w + ((^) ), giving a small 

1/rg contribution. 

A summary of SEst is given in Tab. |T] For e = 0, 
a sweet spot is present for both the (1,1) and (3,3) 
cases. SEst comes from a direct coupling of TLFs to 
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Table I. Influence of CTs on two- and six-electron DQDs. 
Est is shifted, depending on the distance ro between CT and 
DQD. Two sweet spots e = and e = e HB are identified (cf. 
Fig. lk. The hybridization factor (js/t /Us/t) (parameter 
introauced in Eq. |l])) enhances the coherence time for e = 0. 
SEst decreases with ro. Note that for the six-electron DQD 
the decay is much faster at e = tHB- the CTs and the qubit 
couple only in second-order perturbation theory. 



the quadrupole moment of the DQD. The energy shifts 
are on the order of a few gigahertz, corresponding to a 
dephasing time of ns. This time scale is consistent with 
experiments on DQD charge qubits |3"5l 155] . Another 
sweet spot is identified at e = £hb- For two-electron 
DQDs the scaling in ro is identical to e = 0, only lacking 

the hybridization factor (ts/t/Us/t) ■ Ti is improved 
for six-electron DQDs, as the CTs modify Est coupling 
only to the quadrupole moment in second order. 

Robust Single Qubit Gating — We have identified 
two points e = and e = €hb that are well isolated 
from external noise sources. It is possible to manipu- 
late the qubit while staying mainly at these sweet spots. 
Changing the magnitude of Est produces a phase gate: 
U = Ja z , J = f dr Est(t). Est is small at e = 0, 
while at e = £hb, Est = A. A possible gate sweep 
starts from e = 0, tunes the bias rapidly to e = €hb', 
after some waiting time the bias is brought back to e = 
(cf. inset of Fig. [2]). While the manipulation must be 
fast to avoid charge noise, it should still be adiabatic 
with respect to the coupling to excited states (cf . Fig. [TJ . 
The slew rate is limited by the leakage to higher states, 
which is approximated with the transition probability at 
a Landau-Zener crossing of strength r which is crossed 

2 ^ 

with velocity v s i ew [37]: Plz = e . Since the 

tunnel coupling enters Plz quadratically, realistic values 
of t allow very fast manipulations with permitted pulse 
lengths far below nanoseconds. 

We show a fidelity analysis of a 7r-phase gate for a two- 
and six-electron DQD in Fig. [2j The slew rates are fixed 
through Ptz to produce negligible leakage [3S]. We use 
similar densities of the CTs for the two- and six-electron 
DQDs, which are positioned randomly around the DQD 
to generate 1// noise; the coupling to CTs vary the pa- 
rameter A through electrostatic couplings to the DQD 
potential. We exclude a volume around the QD, where 
no CTs are permitted; such nearby TLFs make the DQD 
completely nonfunctional as a qubit. We take the ex- 
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Figure 2. Fidelity analysis for 7r-phase gate for a two/six- 
electron DQD, shown in red/blue. Points are from simula- 
tions involving 1/f noise sources. The fidelities are poor for 
slow manipulation times, which are required by small tunnel 
couplings r; cf. the transition probability at a Landau-Zener 
crossing Plz- Increasing r allows faster qubit manipulations, 
which increases the fidelity. The fidelity of the six-electron 
DQD approaches 1, while it stays much lower for the two- 
electron system. The solid lines are fits using Eq. with 
T2 — 1.5/29.3 ns for the two/six-electron system. The inset 
describes the pulse profile of a 7r-phase gate. Starting from 
e — 0, the DQD is biased to e = ens\ we linearly increase 
e for a time t a i ew . The qubit stays at e = ens for t wa it\ fi- 
nally the qubit is brought back to e = 0, picking up in total 
an odd number of 7r rotations. The overall gate time equals 

2£sZesu ~t~ t-wait- 



eluded volume for two-electron DQDs to be considerably 
larger than for the six-electron system. Fluctuations in 
the tunnel coupling or the pulse profiles are disregarded. 
The sweet spots, especially e = 6hb, offer the advantage 
that Est does not change over a wide range of e. 

The fidelity of the gate, both for the two-electron and 
the six-electron systems (blue/red), is low for small tun- 
nel couplings t. The fidelity increases very quickly with 
r for six-electron DQDs and reaches an ideal value very 
close to 1 . The improvement of the fidelity for the two- 
electron system is much slower. We approximate the 
curves according to Eq. ^ , yielding a coherence time of 
1.5 ns for the two-electron system and 29.3 ns for the six- 
electron case. Steps seen for the two-electron system are 
generated by different waiting times in (2, 0) when con- 
structing a 7r-phase gate; a one-parameter fit to Eq. Q 
cannot completely reproduce these results. 

Conclusion — We propose a fast and robust way to 
manipulate singlet-triplet qubits (STQs) via a high-bias 
phase gate. Contrary to current realizations of phase 
gates, our approach works by going to high bias. The 
qubit couples weakly to CTs; we manipulate rapidly be- 
tween two sweet spots. The "high bias" sweet spot chb 



is not at a specific point in the charge diagram; there is 
a large range of parameters where Est is constant. Note 
that the Rabi rotation gate needed for full qubit control 
is envisioned to occur also at a sweet spot (at e = 0), 
employing magnetic field gradients. It is worth pointing 
out that the proposed high-bias phase gate works also 
as a maximally entangling two-qubit gate for single QD 
qubits p]. 

It would be favorable for our proposal for DQDs to 
have small singlet-triplet energy splitting at e = tuB 
(A, cf. Fig. [IJ, to give comfortable electrical manipula- 
tion times (sub-nanosecond has become accessible [39 J. 
DQDs with A on the order of 30 GHz have been reported 
[T4] . One can decrease the singlet-triplet energy splitting 
further by using favorable dot sizes and external mag- 
netic field parameters. Indeed, we note that a transition 
from a singlet to a triplet ground state is indicated in 
calculations on four-electron QDs |40j . 

A clear prediction of our work is that the many- 
electron QDs, specifically those for which the valence 
electrons occupy the second shell, are uniquely suited to 
protect STQs from charge noise. The mechanism we use, 
resulting from the equality of singlet and triplet charge 
distributions in the second shell, is distinct from pro- 
posals to use screening effects in multi-electron DQDs 
|41l |4"2"] . The manipulation of our six-electron STQs can 
be performed the same way as for the two-electron DQDs, 
including initialization, manipulation, and measurement. 
Additional noise sources, which couple in via the charge 
density, like pure phonon dephasing |431 HI] , are also di- 
rectly suppressed in our approach. We are hopeful that 
the prospect of an order of magnitude improvement in 
gate fidelity will motivate the further experimental ex- 
ploration of the multi-electron regime in quantum dot 
qubits. 
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1 Description of Fidelity Analysis 

We model charge noise acting on double quantum dots (DQD) by a random distribution of charge traps, being 
either filled or empty (cf. Fig. 1). We simulate the time evolution of the DQD numerically using quantum process 
tomography [1]. We generate randomly a distribution of two level fluctuators (TLF), with a broad range of switching 
rates 7. A reasonable probability distribution is P (7) ~ I/7 [2]. The charge distribution is constant during one run 
of the simulation, while the potential fluctuates between successive simulations. This scenario mimics consecutive 
single shot measurements, with a long time between the measurements. 

The coupling strength of DQD and TLF is determined by their distance. As described in the main text, we take 
the shift in the singlet-triplet splitting SEst as the only dynamic variable. For the (2,0) and (4,2) configurations, 
the energy shifts are: 

^^(^■^(^-^^((W) (D 

The excessively occupied QD is positioned in the x-y plane at the coordinate origin, while charge traps are randomly 
occupying the space around the DQD. 

We use material parameters of GaAs. The confining strength Hujq = 3 meV is a common approximation for 
QDs [3]. w c /^o = 0.1 describes moderate external magnetic fields of 0.7 T. The singlet-triplet splitting A is rather 
small, consistent with Dial et al. [4]. All parameters are summarized in Tab. 1. 



The electron distribution can be approximated by the spread of the ground state wave function: ob = 




20 nm. We use 250 TLF with a distance [2.5, 15] as from the coordinate origin for the six-electron system. For 
the two-electron system, we need to exclude a larger volume around the DQD. Otherwise the energy shifts due to 
Eq. (1) destroy the qubit fidelity completely. To generate the same density of TLFs around the DQD, we include 
196 charge traps with a distance [15, 25] as from the origin. 
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Figure 1: Model of charge traps which couple electrostatically to a DQD. The charge traps fluctuate between being 
filled with one electron or being empty. The charge traps are randomly generated for a single run of the simulation. 
They fluctuate slowly; during one gate simulation their electrostatic potential is constant. 



Name 


Value 


energy difference between 
(1,1)/ (2,0) or (4, 2)/ (3, 3) Ua 
singlet-triplet splitting A 

dielectric constant e r 
effective mass m 
confining energy fiwo 
magnetic energy Huj c 


0.5 mcV 
10 GHz 

12.5 
0.067 • TO e 
3 meV 
O.lhujo 



Table 1: Parameters used for the simulation of the STQs. Us and A are chosen to describe the DQD dynamics 
according to Eq. (1) in the main text. The dielectric constant e r and effective mass m corresponds to GaAs: u 
and uj c mimic common confining strengths and magnetic fields of 0.7 T. 
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